library(pracma)
###############################################################
# x=c(9,1,8,2,7,3,6,4,5)
# 
# s=0
# for (k in 1:(length(x)-1)){
#   for (m in (k+1):length(x)){
#     if (x[k]>x[m]) s=s+1
#   }
# }
# print(s)


###############################################################
A=matrix(c(1,2,3,20,4,5,4,39,3,2,1,16),nrow=3,byrow=T)
print(A)
print(rref(A))





